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Concept of spectral soliton coupling is proposed to explain the group-velocity-matching dispersive wave gen- 
eration and the soliton spectral tunneling effect. Soliton eigen state is found corresponding to a spectral phase 
profile under which the soliton phase shifts induced by the dispersion and the nonlinearity are instantaneously 
counterbalanced. A local but non-eigen soliton pulse will shed off energy and form radiations at other wave- 
lengths, during the coupling to the soliton eigen state. Group- velocity-mat hing between the local and generated 
wave is necessary not only as it produces a coupler-like phase profile which supports soliton coupling but also 
because it physically supports a soliton phase matching condition, which means the generated wave has a soliton 
state phase matched to the local exciting pulse. Examples in realistic photonic crystal fiber structures are also 
presented. 
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Dispersive waves (DWs), also known as optical 
Cherenkov radiations (OCRs) 1 , are generated when 
temporal solitons are perturbed by higher-order disper- 
sions [2]. This phenomenon was experimentally verified 
soon after its proposal [3j [4j . DWs play an important 
role in the attractive octave-spanning supercontinuum 
generation (SCG) in fiber structures as they dominate 
the blue-shifted edge of the spectrum while the Raman 
induced soliton self-frequency shift (SSFS) leads to the 
red-shifted edge [5] . The mechanism behind is the phase 
matching between the DWs and the local soliton wave, 
which can be shown through a PM topology [6]. 

Based on the DW generation, soliton spectral tunnel- 
ing (SST) effect was proposed [7] and investigated [8j|9] 
as a soliton spectral switching phenomenon, in which 
the DWs are actually generated in another anomalous 
group velocity dispersion (GVD) region and form a soli- 
ton state. Fundamentally, SST is driven by the Raman 
induced SSFS and requires a potential barrier in the 
GVD 7 . A typical GVD barrier is a normal GVD re- 
gion sandwiched by two anomalous ones, which is usu- 
ally realizable in photonic crystal fibers (PCFs), with 
the waveguide dispersion reducing the material disper- 
sion and forming multiple zero-dispersion wavelengths 
(ZDWs). Moreover, our recent investigation on SST [To] 
showed that along with the PM condition, group velocity 
(GV) matching is another significant premise for decent 
soliton switching with high proportion and broad band- 
width. 

In this letter, we point out that the GV-matching DW 
generation and the SST effect can be understood as a 
spectral soliton coupling from a local state to its eigen 
state. Analogous to a spatial waveguide in which eigen 
modes are supported under a certain phase profile (spa- 
tial refractive index profile), for temporal solitons, eigen 
state is also supported under a spectral phase profile and 



the soliton phase shifts induced by both the dispersion 
and the nonlinearity are counterbalanced. DW genera- 
tion and SST effect are usually evoked by launching a 
local soliton state but not the eigen state, therefore spec- 
tral soliton coupling between the local and eigen states 
can occur, which shed off the pulse energy and form 
waves at other wavelegnths. GV-matching condition is 
essential to tell whether the generated wave is soliton or 
radiation as it defines the spectral phase profile. If the 
generated wave is GV-matching to the local wave, the 
phase profile performs like a coupler structure and the 
generated wave can form a soliton state, while with GV 
mismatched, the phase profile is a leaking structure and 
the generated wave performs as leakage from the local 
wave, namely forming radiation. 

We begin with the demonstration of the soliton 
eigen state. The model behind is the nonlinear wave 
equation in frequency domain (NWEF) [II 12 with 
the dispersion and the cubic nonlinearity included, i.e. 
§- z E + ij3E + i^F[\E\ 2 E] =0 , where E is the electric 
field written in frequency domain, /3 refers to the phys- 
ical phase, 7 is the nonlinear coefficient and F indi- 
cates the Fourier transform. Pulse self-steepening is au- 
tomatically included in such a spectral equation. The 
phase profile is defined with the initial phase /3q and 
the global GV v 9i q eliminated from the physical phase, 
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i.e. /3 eig (cj) = f}{u) - A) - oj • 

coupler-like phase profile is demonstrated to support a 
soliton eigen state which consists of two spectral peaks 
corresponding to the coupler channels. 

If only a local soliton is launched into the above phase 
profile, the soliton spectral coupling is observed as shown 
in Fig. [ljb) . The local wave shed off energy into the 
adjacent channel and form a wave. Meanwhile, the local 



soliton will experience a spectral recoil effect 13 away 



from the coupling. The energy coupled out will also be 
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Fig. 1. Soliton eigen state, soliton coupling and soli- 
ton leaking shown through simulating the NWEF. The 
nonlinear coefficient 7 = 10(km • W) _1 . (a) soliton eigen 
state of a coupler-like phase profile, each soliton peak has 
power 0.4kW; (b) soliton coupling from the local state to 
the eigen state. The local soliton has hyperbolic secant 
shape, the soliton order is set to unity, FWHM = 17.63 
fs; (c) SST effect with extra 20% material Raman effects. 
The Raman spectral response is chosen from Ref. [14] ; 
(d) soliton leaking in a leaking phase profile; (e) and (f) 
XFROG pattern of (c) and (d) at half propagation dis- 
tance, with the corresponding spectra and profils of the 
electric field; Gate function in the XFROG pattern has 
a hyperbolic secant shape with FWHM = 17.63 fs. 



back since they have the GV-matching condition, and 
the local soliton will be pulled back and then evokes 
the coupling over again. The generated wave (used to 
be understood as DW) here is actually a soliton wave 
which travels together with the local soliton, it is fed 
by the local soliton periodically but unfortunately the 
energy is low due to the recoil effect. 

Including material Raman effects, SST effect is ob- 
served as shown in Fig. [TJc). Since the Raman response 
induces SSFS which is always red-shifted and could be 
against the recoil effect, the local soliton can be kept 
at its position and the soliton coupling continuously oc- 



curs until the local wave are fully coupled into the ad- 
jacent channel to form a new soliton wave. Afterwards, 
the new soliton will also experience the Raman induced 
SSFS, which means the coupling backwards are greatly 
suppressed since the soliton will be red-shifted away from 
the coupling. The XFROG pattern of the SST also shows 
that the generated soliton has the GV matching with the 
local one, see Fig. [TJe). 

However, in a leaking structure of the phase profile 
where the generated wave has GV mismatched to the 
launched soliton, SST will not occur, shown in Fig.JjJd). 
In the XFROG pattern, waves are gradually radiated as 
the leakage from the local soliton, leaving a long tail in 
the pattern, see Fig. [TJf). 

In fact, the GV-matching condition helps turning the 
DW PM condition into a soliton PM condition. The DW 
PM condition can be written as: 

(3(uj r ) = /3 S ol,oj s (Ur) (1) 

where f3 so i j0Js (uj) = /3(uj s ) + (u - uj s )v~1 + q s represents 
the non-dispersive soliton phase with spectrum centered 
at uo Sl q s is the soliton wave number and, for fundamental 
solitons, its contribution is minimal [15] . Therefore, the 
above equation can be expanded as: 

/3(u r ) + (U ~ 0J r ) • Vg l = P(UJ S ) + ((J- 0J S )V~1 (2) 

With GV-matching condition, v g ^ r = v g ^ si the above 
equation finally becomes f3 so i^ s (u) = f3 so i j0Jr (uj), imply- 
ing the phase matching within the whole frequency do- 
main between two solitons located at uj s and uj r . 

In practice, PCFs with designed pitch A and hole d 
sizes can achieve flexible dispersion profiles with multi- 
ple controlled ZDWs. For example, in Fig. [2] a solid-core 
index-guiding PCF with a triangular air-hole partten in 
the cladding has form a concave-like dispersion profile 
with 3 ZDWs. Such dispersion profile is actually pro- 
duced by the mode coupling between the core and the 
air-hole cladding around a resonant wavelength. It is no- 
ticed that at long wavelengths, the mode field distribu- 
tion is expanded into the cladding and therefore the non- 
linear coefficient is reduced as the effective mode area is 
increased. 

The phase profile of such a PCF structure is plotted in 
Fig. [3|a), by eliminating the /?o and v 9i o around 1.5/im, 
showing a coupler-like profile. Under such a profile, a lo- 
cal fundamental soliton centered at 1.4/im is launched, 
which is little away from the coupling. Hence, during 
the propagation, the local soliton will first experience 
the Raman induced SSFS, red-shifting towards the cou- 
pling position. Then, strong soliton coupling occurs with 
most of the energy transferred into the adjacent channel, 
centered at 1.85/im. Actually, the choose of the reference- 
like and v 9i o is directly linked to the PM condition. A 
symmetry coupler-like phase profile, with the same peak 
height, is corresponding to the PM between the local 
and nonlocal solitons, under which the soliton coupling is 
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Fig. 2. Solid-core index guiding PCF's dispersion and 
effective nonlinear coefficient 7 versus wavelength, esti- 
mated by means of Comsol software. Insets are mode 
field distributions at different wavelengths. 

strongest. The XFROG pattern in Fig.[3jc) again proves 
the GV matching between the two solitons. The propor- 
tion of the soliton energy transfer is around 60%, lower 
than the ideal SST because the nonlinear coefficient 7 
has a reduction in long wavelengths. 

With a different pitch size, the phase profile turns to 
a leaking structure. Therefore, the soliton leaking oc- 
curs instead of coupling, see Fig. [3^b). Although from 
the spectral evolution, the nonlocal wave generated are 
still decent, they are actually not soliton waves but just 
leakage-like DWs, proved by the XFROG pattern in 

Fig. id). 

As a conclusion, we propose the concept of soliton cou- 
pling from a local state to the eigen state, which helps us 
properly explain the DW generation and SST effect. The 
basis of the soliton coupling is the soliton eigen state un- 
der a spectral phase profile, in which the soliton phase 
shifts induced by both the dispersion and the nonlin- 
earity are counterbalanced. Under a coupler-like phase 
profile, soliton eigen state is demonstrated to have a two- 
peak profile, and a local but non-eigen soliton is there- 
fore expected to excite the soliton coupling. Sepctral re- 
coil effect is observed during the coupling, which always 
pushes the pulse away from the coupling position and 
suppresses the coupling. With material Raman effects, 
red-shifted SSFS is produced and could be against the 
recoil effect, which therefore helps keep the pulse at its 
position and causes continuously soliton coupling, known 
as SST effect. The GV-matching condition is crucial not 
only as it produces such a coupler-like phase profile, but 
also because it physically turns a DW PM condition into 
a soliton PM condition. With GV mismatched between 
the local pulse and the generated wave, the phase profile 
is supposed to be a leaking structure, under which the 
generated wave performs as leakage from the local pulse. 
Such soliton coupling and soliton leaking effects in real- 
istic PCF structures are presented, in which SST with a 
spectral span over 300 nm and 60% energy transferring 
is numerically demonstrated. 
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Fig. 3. NWEF simulations on the spectral evolutions 
of a 25 fs (FWHM) local fundamental soliton centered 
at 1.4/im under the coupler-like and leaking phase pro- 
files: (a) coupler-like structure with A = 2.255/im and 
(b) leaking struct uere with A = 2.235/im. The input, 
output spectra of both cases and broadest spectrum dur- 
ing the SST (position marked by white dashed line) are 
shown. The input pulse has a peak power 380W, the 
material Kerr nonlinearity is 2.6 x 10 _20 m 2 /W, Raman 
fraction is 24.5%. (c) and (d) XFROG pattern of (a) and 
(b) at a propagation length of 20m. The Gate function 
in XFROG pattern also has a FWHM = 25 fs. 
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